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^ ' Abstract 

We discuss some challenging open problems in the geometric con- 
^^ I trol theory and sub-Riemannian geometry. 
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It is getting harder to prove theorems and easier to force other people 
to prove them when you are sixty. Some colleagues asked me to describe 
interesting open problems in geometric control and sub-Riemannian geom- 
etry. Here I list few really challenging problems; some of them are open 
for a long time and were publicly or privately stated by well-known experts: 
^ '. J.-M. Coron, I. Kupka, R. Montgomery, B. Shapiro, H. Sussmann, and oth- 
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CN . I- Singularities of time-optimal trajectories. 



ers. 



Let f,g he a pair of smooth (i.e. C°°) vector fields on a n-dimensional 



cn ' manifold M. We study time-optimal trajectories for the system 

Q = fici)+ug{q), |n| < 1, g G M, 



with fixed endpoint. Admissible controls are just measurable functions and 
admissible trajectories are Lipschitz curves in M. We can expect more regu- 
larity from time-optimal trajectories imposing reasonable conditions on the 
pair of vector fields. 

A. (/, g) is a generic pair of vector fields. Optimal trajectories cannot 
be all smooth; are they piecewise smooth? This is true for n = 2. More 
precisely, if dimM = 2, then any point of M has a neighborhood such 
that any contained in the neighborhood time-optimal trajectory is piecewise 
smooth with atmost 1 switching point (see [321 HI]). According to the control 
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theory terminology, a switching point of an admissible trajectory is a point 
where the trajectory is not smooth. 

The question is open for n = 3. What is known? Let sw(q) be mini- 
mal among numbers k such that any contained in a sufficiently small neigh- 
borhood of g G M time-optimal trajectory has no more than k switching 
points. We set sw(g) = oo if any neighborhood of q contains a time-optimal 
trajectory with an infinite number of switching points. It is known that 
sw(g) = 2 for any q out of a 2-dimensional Whitney stratified subset of the 
3- dimensional manifold M (see [331 I2H]) and sw < 4 for any q out of a 1- 
dimensional Whitney stratified subset of M (see [7]). Some further results 
in this direction can be found in [2Ij. We do not know if sw(g) < oo for 
any q G M. We also do not know if a weaker property, the finiteness of 
the number of switching points for any individual time-optimal trajectory is 
valid. 

Higher dimensions. There is a common opinion that starting from some 
(not very big) dimension, time-optimal trajectories with accumulating switch- 
ing points cannot be eliminated by a C°°-small perturbation of the system 
and thus survive any genericity conditions. However, to my knowledge, this 
opinion was never supported by a proof. There are very interesting examples 
of extremals with accumulating switching points whose structure survives 
small perturbations (see [211 E3) but nobody knows if these extremals are 
optimal. 

B. /, g are real analytic vector Gelds. Let M be a real analytic manifold 
and /, g analytic vector fields, not necessary generic. Here we cannot expect 
any regularity of an arbitrary time-optimal trajectory. Indeed, it is possible, 
even for linear systems, that all admissible trajectories are time-optimal. We 
can however expect that among all time-optimal trajectories connecting the 
same endpoints there is at least one not so bad. 

If n = 2, then any two points connected by a time-optimal trajectory can 
be connected by a time-optimal trajectory with a finite number of switching 
points (see [3S1 ES])- This is not true for n > 3. Indeed, classical Fuller 
example with accumulating switching points [17] can be easily reformulated 
as a 3-dimensional time-optimal problem. Main open question here is as 
follows: Given two points connected by a time-optimal trajectory, can we 
connect them by a time-optimal trajectory with no more than a countable 
number of switching points? 

What is known? The points can be connected by a time-optimal trajec- 



tory whose set of switching points is nowhere dense [M], is not a Cantor set 
(can be derived from pQ), and satisfies some additional restrictions |31|. We 
do not know if we can avoid a positive measure set of switching points. 

All mentioned open questions are not easy to answer. In my opinion, the 
most interesting is one on generic 3-dimensional systems. 

II. Cutting the corners in sub-Riemannian spaces. 

Unlikely the just discussed problems, optimal paths in sub-Riemannian 
geometry are usually smooth. However we do not know if they are always 
smooth. A natural open question here is as follows. Let 7i : [0, 1] — )■ M be 
two smooth admissible paths of a sub-Riemannian structure on the manifold 
M, 7o(0) = 7i(0) = go, 7o(0)A7i(0) 7^ 0. Does there exist an admissible path 
connecting 7o(l) with 71 (1) that is strictly shorter than the concatenation of 
the curves 70 and 71? 

Admissible paths are integral curves of a bracket generating vector dis- 
tribution A C TM. It is easy to show that positive answer to the question 
for rank 2 distributions implies positive answer in the general case. Let 
A = span{fQ, /i}, where /o, /i are smooth vector fields on M, dimM = n. 
We set 

nk{q) = span{[fi^, [■ ■ ■ , /.J ■ ■ ■]{q) : ij G {0, 1}, j < k} , 

m = min{A; : nk{qo) = n}. li m < A, then the answer to our question is 
positive: it is proved in [22]. Moreover, an example studied in [21] supports 
the conjecture that the answer is perhaps positive for ?7i = 5, n < 4 as well. 
Any improvement of the estimates for m and n would be very interesting. 
We still know very little about sub-Riemannian structures with big m and it 
may happen that the answer is negative for some m and n. 

III. "Morse— Sard theorem" for the endpoint maps. 

We continue to consider admissible paths of a sub-Riemannian structure 
on M. Given go ^ M, the space of starting at go admissible paths equipped 
with the if ^-topology forms a smooth Hilbert manifold. The endpoint map is 
a smooth map from this Hilbert manifold to M ; it sends a path 7 : [0, 1] — )■ M 
to the point 7(1). Critical points of the endpoint map are called singular 
curves of the distribution. 

The Morse-Sard theorem for a smooth map defined on a finite dimen- 
sional manifold states that the set of critical values of the map has zero 



measure. It is not true in the infinite dimensional case: there are smooth 
surjective maps without regular points from any infinite dimensional Banach 
space to M^ (see j9]). 

The endpoint maps have plenty of regular points but we do not know if 
they always have regular values. This is an interesting open question. We 
can reformulate the question as follows: is it possible that starting from gg 
singular curves fill the whole manifold Ml 

Optimal (i. e. length minimizing) singular curves are better controlled; 
we know that starting from go optimal singular curves fill a nowhere dense 
subset of M (see [2j). An important open question: can they fill a positive 
measure subset of Ml 

IV. Unfolding the sub-Riemannian distance. 

The problem concerns singularities of the distance function for generic 
sub-Riemannian structures. Let q^ & M and Sq^ : M ^ [0, +cxd) be the sub- 
Riemannian distance from the point go- Sufficiently small balls ^^^([O,^]) 
are compact. Let g G M be a point from such a ball. Then g is connected 
with go by an optimal path. If this optimal path is not a singular curve and 
any point from a neighborhood of g is connected with go by a unique optimal 
path, then Sq^ is smooth at g. 

The points connected with go by more than one optimal path form the 
cut locus. The function Sq^ is not smooth in the points of the cut locus and 
it is not smooth at the points connected with go by optimal singular curves 
but these two types of singularities are very different. 

If all connecting go and g optimal paths are not singular curves, then 
the singularity of Sq,^ at g is similar to singularities of Riemannian distances 
and, more generally, to singularities of the optimal costs of regular variational 
problems. The function S^q is semiconcave P^ and typical singularities in low 
dimensions are well-described by the theory of Lagrangian and Legendrian 
singularities [Sj Ch. 3] developed by V. Arnold and his school. 

On the other hand, if go is connected with g by an optimal singular curve, 
then SqQ is not even locally Lipschitz at go (see ^ Ch. 10]); moreover, classical 
singularities theory does not work and the structure of typical singularities is 
totally unknown. There are few studied models [H [27] but they are too sym- 
metric to be typical and the structure of their singularities is easily destroyed 
by small perturbations. 

Let us consider, in particular, the Martinet distribution that is a rank 
2 distribution in M? in a neighborhood of a point go such that n2(go) = 



2, n^igo) = 3 (see II. for the definition of nj(go))- The points g in a neigh- 
borhood of go where n2{q) = 2 form a smooth 2-dimensional submanifold 
A^ C M, the Martinet surface. Moreover, the distribution A is transversal 
to A^ and A^ fl TgN, q & N is a. fine distribution on N. Integral curves of 
this line distribution are singular curves whose small segments are optimal. 
There are no other singular curves for such a distribution. 

Example. Let /i = ^, /a = ^ + xj-£^, then A = span{fij2} is a 
Martinet distribution and the Martinet surface is a coordinate plane defined 
by the equation xi = 0. The fields /i, /2 form an orthonormal frame of the so 
called 'flat' sub-Riemannian metric on the Martinet distribution. Let go = 0, 
singularities of 5*0 are well-known (see [1]). The cut locus has the form: 
{x G M^ : Xi = 0, X2 7^ 0}, the Martinet surface with the removed singular 
curve through go. The singular locus of a sphere S~^^{e) is a simple closed 
curve and its complement (the smooth part of the sphere) is diffeomorphic 
to the disjoint union of two discs. 

The 'flat' metric is rather symmetric, in particular, it respects the orthog- 
onal reflection of M^ with respect to the Martinet plane. Simple topological 
arguments show that for generic metric with a broken symmetry, the smooth 
part of a sphere is connected and is not contractible. The singular locus of 
the sphere should be cut at the points where the sphere intersects the optimal 
singular curve but the shape of the sphere near these points is unknown. 

An important open question is to flnd a C^-classification of the germs of 
spheres at the points of optimal singular curves for generic metrics. Here we 
say that two germs are C^-equivalent if one can be transformed into another 
by a germ of C^-diffeomorphism of M^. 

The next step is the Engel distribution, i. e. a rank 2 distribution in M^ 
such that n^lq) = 3, n^lq) = 4. There is exactly one singular curve through 
any point and small segments of singular curves are optimal. We repeat our 
question for this case; the spheres are now 3-dimensional hypersurfaces in 
R\ 

Generic germ of a rank 2 distribution in R" possesses a (n— 4)-dimensional 
family of singular curves through go G M" (see, for instance, [25]). Take a 
generic curve from this family; its small segments are optimal. Take a point 
g where the selected singular curve intersects the sphere S~^^{e). The points 
of singular curves from the family in a small neighborhood of g in our sphere 
form a smooth {n — 4)-dimensional submanifold S C M". The intersection 
of the sphere with a transversal to S smooth 4-dimensional submanifold 
should have a shape similar to the germ of the sphere in the Engel case. A 



neighborhood of q in the sphere is fibered by such intersections. Hence the 
solution of the problem in the 4-diniensional Engel case is a very important 
step in the unfolding of the sphere for any n > 4. 

The desired classification seems to be complicated. There is a 2-dimensional 
modification of the problem that, in my opinion, already contains the main 
difficulty. Ones resolved, it will reduce the study of higher dimensional prob- 
lems to the conventional singularities theory techniques. Consider the germ 
at go G M^ of a pair of smooth vector fields /o, /i such that fo{qo) A /i(go) = 
0, n2{qo) = 1, nsi^qo) = 2. The almost Riemannian distance Sqg{q) is the 
optimal time to get q from go by an admissible trajectory of the system 

q = Uofoiq)+Uifi{q), u\^u\ = \. 

The question is to find a C^-classification of the germs of distance functions 
S'go for generic pairs of vector fields /o, /i among the pairs that satisfy con- 
ditions 71-2 (go) = 1, ^^3(^0) = 2. See [10] for some partial results. 

V. Symmetries of vector distributions. 

A symmetry of a distribution A C TM is a diffeomorphism $ : M — )■ 
M such that <I'*A = A. The differential geometry appeals to search most 
symmetric objects in the class, those with a maximal symmetry group. The 
singularities theory, on the contrary, encourages the study of less symmetric 
generic objects. Both paradigms have their reasons and complement each 
other. Anyway, a fundamental problem is to characterize objects whose 
symmetry groups are finite-dimensional Lie groups. 

Our objects are vector distributions. Any symmetry transfers singular 
curves of the distribuiton in singular curves and these curves often play a 
key role in the calculation of symmetry groups (see [2D1 US])- We say that a 
distribution is singular transitive if any two points of M can be connected by 
a concatination of singular curves. A natural open question is as follows: Is 
it true that singular transitivity of the distribution implies that its symmetry 
group is a finite dimensional Lie group? 

All known examples support the positive answer to this question. More- 
over, the group of symmetries is infinite dimensional for many popular classes 
of not singular transitive distribution: codimension 1 distributions, involu- 
tive distributions, Goursat-Darboux distributions. We can even expect that 
any not singular transitive rank 2 distribution has an infinite dimensional 
symmetry group. Some results of [15] seem to be rather close to this state- 
ment. 



VI. Closed curves with a nondegenerate Frenet frame. 

Let 7 : S*^ — 7- M" be a smooth closed curve in MJ^. We say that 7 is degen- 
erate at t G S*^ if 7(t) A ■ ■ ■ A 7^"''(t) = 0. Degeneracy points are the points 
where velocity or curvature of the curve vanishes if n = 2, where velocity or 
curvature or torsion vanishes if n = 3 e. t. c. The curve is nondegenerate if it 
has no degeneracy points. Any nondegenerate curve admits the orthonormal 
Frenet frame E{t) = (ei(t), . . . , e„(t)) , t E S^, that is a smooth closed curve 
in the orthogonal group 0(n). 

Now let n = 3 and 7 be a plane convex curve, 7(t) G M^ C M.^, Vt G 5*^. 
Then any small perturbation of 7 as a spatial curve is degenerate in some 
points. On the other hand, an appropriate small perturbation of a plane 
convex curve run twice (say, of the curve t 1— )■ 7(2t), t G S^) makes it a 
nondegenerate curve in M.^. Everyone can get evidence of that playing with 
a cord on the desk. This is also a mathematical fact proved in [161 123] . 

Frenet frame of the plane convex curve treated as a spatial curve is a 
one-parametric subgroup SO (2) C 0(3), a shortest closed geodesic in 0(3) 
equipped with a standard bi-invariant metric. It is proved in [23] that the 
length of the Frenet frame of any regular curve in M'^ is greater than the 
double length of SO (2). 

Come back to an arbitrary n. Let fi{n) be minimal m such that run m 
times convex plane curves have regular small perturbations in M". We know 
that /x(2) = 1, yu(3) = 2. An important open problem is to find fi{n) for 
n > 3 and to check if the length of the Frenet frame of any regular curve in 
M" is greater than the length of SO (2) C 0(n) multiplied by fi{n). 

Let me explain why this problem is a challenge for the optimal control 
theory and why its study may bring important new tools to the theory. The 
Frenet structural equations for a regular curve 7 in M" have a form: 

7 = ei, Ci = Ukit)ei+i - Mi_i(t)ei_i, i = l,...,n-l, (1) 

where Uq = Un = 0, Uiit) > 0, i = 1, . . . , n — 1, t G S^. 

In other words, regular curves together with there Frenet frames are pe- 
riodic admissible trajectories of the control system (1) with positive control 

parameters mi, . . . , «„. The length of the Frenet frame on the segment [0, ti] 

ti 1 

is / {u\{t) + ■ ■ ■ + M^_i(t)) ^ dt. We are looking for a periodic trajectory with 


shortest Frenet frame. 



A shortest frame is unlikely to exists since control parameters belong to 
an open cone. It is reasonable to expect that minimizing sequences converge 
to a solution of (1) with U2(t) = ■ ■ ■ = Un-i{t) = 0, while ui{t) stays positive 
to guarantee the periodicity of 7. In other words, the infimum is most likely 
realized by a plane convex curve run several times. Obviously, the length of 
the Frenet frame does not depend on the shape of the convex curve. 

So we have to take the m times run circle: Ui{t) = 1, U2(t) = ■ ■ ■ = 
Un-i(t) =0, < t < 2nm, and try to find small positive perturbations of 
control parameters in such a way that the perturbed curve stays periodic. 
Then /i(n) is minimal among m for which such a perturbation does exist. 
Unfortunately, we cannot use typical in geometric control sophisticated two- 
side variations that produce iterated Lie brackets: only one-side variations 
are available. I think, it is a very good model to understand high order effects 
of time-distributed one-side variations. 

The study of the 3-dimensional case by Milnor in |23] was not variational; 
it was a nice application of the integral geometry. However, the integral 
geometry method is less efficient in higher dimensions (see [2S] for some 
partial results). 

VII. Controllability of the Navier— Stokes equations controlled 
by a localized degenerate forcing. 

We consider the Navier-Stokes equation of the incompressible fluid: 

du 

— + {u, V)u - uAu + Vp = r]{t, x), divu = 0, (2) 

with periodic boundary conditions: x G T^/27rZ'^, d = 2,3. Here u{t, x) G M"^ 
is the velocity of the fluid at the point x and moment t; i/ is a positive 
constant (viscosity), p is the pressure and rj external force. 

We treat (2) as an evolution equation in the space of divergence free 
vector fields on the torus T*^ controlled by the force. In other words, u{t, ■) 
is the state of our infinite dimensional control system and r^ is a control. 
These notations are against the control theory tradition where u is always 
control but we do not want to violate absolutely standard notations of the 
mathematical fiuid dynamics. By the way, symbol u for the control was 
introduced by Pontryagin as the first letter of the Russian word "upravlenie" 
that means control. 

The state space is \^ = |u G H^{T'^,M.'^) : divu = O}, control parameters 
ri{t, ■) belong to a subspace E (Z V. We say that the system is approximately 
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controllable (controllable in finite dimensional projections) in any time if for 
any Uo,Ui & V, ti > and any e > (any finite dimensional subspace F C 
V) there exists a bounded control r/, 'r]{t, ■) & E, < t < ti, and a solution 
u of (2) such that u(0, ■) = uq, \\u{ti, ■) — mi||l2 < ^ (-Pf(m(^i, ■) — ^i) = 0, 
where Pp is the L2-orthogonal projector on F) . 

Of course, controllability properties depend on the choice of the space of 
control parameters E. It is known that the systems is controllable in both 
senses by a localized forcing when E = {u & V : supp u C V} and V is an 
arbitrary open subset of T''. Moreover, such E provides a much stronger 
exact controllability (see pj fTi| [T8| fT9]). 

On the other hand, the system is approximately controllable and con- 
trollable in finite-dimensional projections by a degenerate forcing (or forc- 
ing with a localized spectrum) when ii^ is a finite dimensional space of low 
frequency trigonometric polynomials (see O O [291 EO]). This kind of con- 
trollability illustrates a mechanism of the energy propagation from low to 
higher frequencies that is a necessary step in the long way towards a reliable 
mathematical model for the well-developed turbulence. 

It is important that the control parameters space E does not depend 
on the viscosity u. Moreover, if d = 2, then the described controllability 
properties are valid also for the Euler equation (i. e. for u = 0); the Cauchy 
problem for the Euler equation is well-posed in this case. 

Now an important open question: is the system approximately control- 
lable and (or) controllable in the finite dimensional projections by a local- 
ized degenerate forcing when ii^ is a finite dimensional subspace of the space 
{u E V : suppw C V}7 The question is about existence and effective con- 
struction of such a space E that does not depend on the viscosity u. 

The independence on u is important for eventual applications to the well- 
developed turbulence that concerns the case of very small u (or very big 
Reynold number). Of course, similar problems for other boundary conditions 
and other functional spaces are also very interesting. 

We have arrived to a sacral number of seven problems and can relax a 
little bit. To conclude, I would like to discuss one more problem; it is less 
precise than already stated questions but, to my taste, is nice and fascinating. 
The problem concerns contact 3-dimensional manifolds and is inspired by the 
Ricci fiow story. 



oo. Diffusion along the Reeb field. 

I recall that a contact structure on a 3-diniensional manifold M is a 
rank 2 distribution A C M such that ^2(5) = 3, Vg G M. According to a 
classical Martinet theorem, any orientable 3-dimensional manifold admits a 
contact structure. I am going to introduce some dynamics on the space of 
sub-Riemannian metrics on a fixed compact contact manifold (M, A). 

First, to any sub-Riemannian metric on (M, A) we associate a transversal 
to A Reeb vector field e on M . In what follows, we assume that A is oriented; 
otherwise e(g) is defined up-to a sign but further considerations are easily 
extended to this case. Let a; be a nonvanishing differential 1-form on M that 
annihilates A. The condition 77,2(9) = 3 is equivalent to the inequality Uq A 
dqUJ 7^ 0. The form u is defined up-to the multiplication by a nonvanishing 
function; the sign of the 3-form ojq A dqU does not depend on the choice 
of OJ and defines an orientation on M. We have: dqu], 7^ 0; moreover, 

dq{auj)\^ = a{q)dqUj\^ for any smooth function a of M. 

Given a sub-Riemannian metric on A, there exists a unique annihilating 
A form oj such that the 2- form dqijj\. coincides with the area form on A„ 
defined by the inner product and the orientation. The kernel of dqU is a 
1-dimensional subspace of TqM transversal to A^, and e(g) is an element of 
this kernel normalized by the condition {uq,e{q)) = 1. 

In other words, the Reeb vector field is defined by the conditions: 
ieOJ = 1, iedu = 0. Hence L^OJ = 0, where Lg is the Lie derivative along 
e, and the generated by e flow on M preserves u. In general, this flow does 
not preserves the sub-Riemannian metric. We may try to classify contact 
structures by selecting best possible sub-Riemannian metrics on them. 

Assume that there exists a metric preserved by the flow generated by the 
Reeb vector field. Take a standard extension of the sub-Riemannian metric 
to a Riemannian metric on M: simply say that e is orthogonal to A and 
has length 1. The generated by e flow preserves this Riemannian metric as 
well. So our compact Riemannian space admits a one-parametric group of 
isometries without equilibria. Hence M is a Seifert bundle. Do not care if 
you do not remember what is Seifert bundle: it is sufficient to know that 
they are classified as well as invariant contact structures on them. 

The invariant with respect to the Reeb field sub-Riemannian metric gives 
a lot of information about the manifold. Let q G M; our sub-Riemannian 
metric induces a structure of Riemannian surface on a neighborhood of q 
factorized by the trajectories of the local fiow generated by the restriction of e 
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to the neighborhood. Let K,{q) be the Gaussian curvature of this Riemannian 
surface at the point q; then /t is a well-defined smooth function on M, a 
differential invariant of the sub-Riemannian metric. Moreover, k is a first 
integral of the flow generated by the Reeb field e. If /t = 0, then universal 
covering of the sub-Riemannian manifold is isometric to the Heisenberg group 
endowed with the standard left-invariant metric. If /t is a negative (positive) 
constant, then universal covering of the sub-Riemannian manifold is isometric 
to the universal covering of the group SL(2) (group SU(2)) equipped with a 
left-invariant sub-Riemannian metric induced by the Killing form. 

Assume that function k, is not a constant and c G M is it regular value. 
Then k,^^{c) is a compact 2-dimensional submanifold of M; we treat it as 
a 2-dimensional Riemannian submanifold of the Riemannian manifold M 
equipped with the standard extension of the sub-Riemannian structure. It 
is easy to see that /t~^(c) is isometric to a flat torus. Indeed, T{k,~^{c)) 
contains the field e|K-i(c) and is transversal to A; the field e|K-i(c) and the 
unit length field from the line distribution T{k,^^{c)) fl A commute and form 
an orthonormal frame. 

So preserved by the Reeb field sub-Riemannian metrics have plenty of nice 
properties. Unfortunately, not any compact contact manifold admits such a 
metric because not any compact 3-dimensional manifold admits a structure 
of Seifert bundle. I am going to discuss a natural procedure that may lead 
to a generalized version of such a metric with reasonable singularities. 

It is more convenient to work in the cotangent bundle than in the tan- 
gent one. A sub-Riemannian metric is an inner product on A C TM; let 
us consider the dual inner product on A* = T*M/A-^, where A-*- is the 
annihilator of A. This is a family of positive definite quadratic forms on 
A* = T*M/A^, q G M, or, in other words, a family of nonnegative quadratic 
forms hq on T*M such that ker hg = A^. The function 

h : T*M -^ R, where h{0 = hg{^), V^ G T*M, q e M, 

is the Hamiltonian of the sub-Riemannian metric. Hamiltonian vector field 
on T*M associated to h generates the sub-Riemannian geodesic flow. 

The Hamiltonian h determines both the vector distribution and the inner 
product. We denote by Uh '■ T*M — )■ R the Hamiltonian lift of the Reeb field 

e, 

UhiO = {^,eiq)), y^eT^M, q e M, 

and by Ul : T*M — )■ T*M, t G M, the Hamiltonian flow generated by the 
Hamiltonian field associated to u^. The flow U^ is a lift to the cotangent 
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bundle of the flow on M generated by e. Let P^ : M — )■ M be such a flow, 
^^ = e o P^{q), Po{q) = g, q e M; then U^ = Pl^. The flow Pt preserves 
the sub-Riemannian metric if and only if the flow Ul preserves h; in other 
words, if and only if {uh, h} = 0, where {■, ■} is the Poisson bracket. Note 
that h\rp,j^ is a quadratic form and m/i|j.*^ is a linear form, Vg G M; hence 

{uh, h}\^^j^ is a quadratic form. 

Recall that the flow U^ = Pl^ preserves the 1-form u, and w is a non- 
vanishing section of the line distribution A-*-. Hence A^ is contained in the 
kernel of the quadratic forms h o Uj^l and {uh, {■ ■ ■ {uh, ^} ■ ■ ■ }L, a^ ~ 

i 
JfL=0 (^°^^) I T,*M- 

We are now ready to introduce the promised dynamics on the space of 
sub-Riemannian metrics on A, where metrics are represented by their Hamil- 
tonians. Let 5 be a positive smooth function on M. A discrete time dynam- 
ical system transforms a Hamiltonian /i„ into the Hamiltonian 

e 

hn+i = — hnoUl^dt, n = 0,1,2,..., 

a partial average of /i„ with respect to the flow W^. 

The Hamiltonian /i„+i is equal to /i„ if and only if {uh„, h^} = 0. Indeed, 
let (■, ■)q be an inner product in A* and if* : A* — )> A* the symmetric 
operator associated to the quadratic form /i„ °W^„|a* ^J ^^^^ inner product: 

hn o ^hS') ~ {^n'y ')q- Recall that the flow Ul^ is generated by the Reeb 
field of hn, hence the area form on A* defined by hnoUj^l does not depend 

on t; in other words, det H!^ = const. The equation det H = const defines 
a strongly convex hyperboloid in the 3-dimensional cone of positive definite 
symmetric operators on the plane, and ii^ is a curve in such a hyperboloid; 

e 

hence -^ J H^dt = H^ if and only if iJ* = H^, i. e. /i„ o U^^ = hn- 

—e 

If the sequence /i„ converges, then its limit is the Hamiltonian of a sub- 
Riemannian metric on A preserved by the Reeb field. Otherwise we may 
modify the sequence and take scaled averages: 

h„j-i = Cn / h„ o Uh dt. 



I'n+l — ^n I if-n ^ t^hn 
-en 
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There is a good chance to arrive to a nonzero hmiting Hamiltonian h^o by 
a clever choice of the sequences of positive functions £:„,c„. Then h^\^^^ 

is a nonnegative quadratic form and A;^ C ker/iooL,»^ for any q E M. It 

may happen however that rank ( /iooL*^ ) < 2 for some q E M and h is not 

the Hamiltonian of a contact sub-Riemannian metric; we can treat it as a 
generalized version of such a metric. 

A continuous time analogue of the introduced dynamics is a "heat along 
the Reeb field" equation 

— = c{uh,{uh,h}} 

in the space of sub-Riemannian metrics on the given contact distribution. 
It is easy to show that the equality {uh,{uh,h}} = implies {uh,h} = 
and stationary solutions of this equation are exactly the metrics preserved 
by the Reeb fields. I conclude with an explicit expression for this nice and 
mysterious evolution equation in the appropriate frame. 

All contact distributions are locally equivalent according to the Darboux 
theorem. Let /i, /2 be a basis of the contact distribution A such that /i, /2 
generate a Heisenberg Lie algebra: [/iJ/1,/2]] = [/2,[/2,/i]] = 0. We set 
Vi{^) = (^, /i(g)), ^ E T*M, qE M, the Hamiltonian lift of the field /i, i = 
1,2; then 

{t;i,{t;i,t;2}} = {t;2,{^^2,t^i}} = 0. (3) 

Hamiltonian of any sub-Riemannian metric on A has a form: 

h = auvl + 2ai2t;it;2 + 0,22^1, (4) 

where aij are smooth functions on the domain in M where f\, /2 form a basis 
of A, and the quadratic form defined by the matrix A{q) = ( "\ (^j a"(a) ) ^^ 
positive definite for any q from this domain. Let S = det A; this is a function 
on M and we treat it as a constant on the fibers function on T*M. A key 
for us function u^ depends only on S and has a form: 

-Uh = 5{Vi,V2} + t^l{t^2, 5} + V2{5, Vi}. (5) 

The relations (3)-(5) give an explicit expression for the equation ^ = 
c{uh,{uh,h}} as a system of third order partial differential equations for 
the functions a^. 
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